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Abstract 

Entangled quantum states are an important element of quantum information techniques. We 
determine the requirements for states of quadrupolar nuclei with spins >l/2 to be entangled. It 
was shown that entanglement is achieved at low temperature by applying a magnetic field to a 
quadrupolar nuclei possess quadrupole moments, which interacts with the electricfield gradient 
produced by the charge distribution in their surroundings. 
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I. INTRODUCTION 

Quantum entanglement [l s| , the most characteristic feature of quantum mechanics, 
is one of the central concepts in quantum information theory and is the feature that dis- 
tinguishes it most significantly from the classical theory. Entanglement is now viewed as 
a physical resource, which provides a means to perform quantum computation and quan- 
tum communication [5] . It should be emphasized that property of entanglement can be 



considered regardless of the nature of qubits 



a. 



One of the most intensively investigated systems is clusters of coupling nuclear spins 



EL 



which received considerable attention as a platform for the practical implementation of 
a quantum computer (QC) by using nuclear magnetic resonance (NMR) technique 
The strength of the coupling (such as dipole-dipole, scalar or exchange interactions) between 
different spins and interaction of the spins with environment determine the decoherence time. 
Short decoherence times limit possible calculation times of the QC [18j. Later it has been 
proposed to used quadrupole nuclei, thus eliminating the requirement of interaction between 



spins and decreasing interaction of qubits with environment 



19 



-|25| 



It was shown that a spin | system is equivalent to a system of two magnetization vectors 
and can be represented using the Pauli spin matrices 2x2 [28;!]. It means that a 
single spin 3/2 is isomorphic to a system consists of two dipolar coupling spins 1/2, which 



can be considered as qubits. 



with quadrupole splitting [21 



'his was experimentally confirmed by using NQR technique 



22| . The feasibility of quantum computing based on a pure 



without external magnetic fields) NQR technique was theoretically investigated in detail 



in 



23]. Using the resonance excitation technique and the level-crossing method, it was 



proposed three quantum logic gates: a controlled NOT, SWAP and NOT 2 . Thus the 
method to synthesize qubits from a set of the spin states of a single particle with spin higher 
1/2 has been developed. It is logical to raise the question of entanglement of these qubits. 

Our present purpose is to investigate entanglement between the quantum states of 
quadrupole nuclei. We consider a nucleus with spin | being in an internal electric field 
gradient (EFG) and an external magnetic field when the quadrupole interaction energy is of 
the order of the magnitude or even greater than the Zeeman one. Results of computer simu- 
lations of entanglement dynamics are presented for real spin systems with a non-equadistant 
energy spectrum. 
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II. NQR IN MAGNETIC FIELD 



Let us consider a quadrupole nucleus with spin / (/ > 1/2) in the thermodynamic 
equilibrium with the density matrix 

Here ks is the Boltzmann constant, T is the lattice temperature, Z = Tr jexp ^— j j is 
the partition function. In the general case, the Hamiltonian "H can consist of the Zeeman 
("Ha/) and the quadrupole ("Hq) parts, 

H = H m + Hq. (2) 

The Zeeman interaction between the applied magnetic field, Hq and nuclear spin is the 
external factor for a crystal and the direction of this field is chosen as the z-axis of the 
laboratory frame, H = H Q z, where H is the strength of the external magnetic field. The 
part of the Hamiltonian describing this interaction is 

U M = -lH Q I„ (3) 

where 7 is the gyromagnetic ratio of the nucleus with spin /, I z is the projection of the 
individual spin angular momentum operators / on the z- axis. 

Quadrupole coupling exists between a non-spherical nuclear charge distribution and an 
electric field gradient (EFG) generated by the charge distribution in their surroundings. It 
is possible to reduce the EFG symmetric tensor to a diagonal form by finding the principal 
axes frame (PAF) with the Z- and X-axises directed along the maximum and minimum of 
EFG, respectively, \V ZZ \ > \V YY \ > \V XX \, where % = (£ = X, Y, Z) and V is the 
potential of the electric field. In the laboratory frame the quadrupolar Hamiltonian can be 
presented in the following form 



eQq. 



4/(27-1) 



H Q = 77%T^ U V><P) 3/ 2 2 -/V?(4-/ 2 ) U+(9,v), (4) 
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where 



U(9 ,(p) =e -i<Plz e -iOIy e i<pIz^ ( 5 ) 

eQqzz is the quadrupole coupling constant of EFG, I± are the raising and lowering operators 
of the spin, and 6 and ip refer to the polar and azimuthal angles determining the orientation 



of the laboratory frame z-axis in the PAF coordinate system. The asymmetry parameter r/ 
is defined as 

Vyy — V X X /„n 
V = 77 , (6) 

Vzz 

which may vary between and 1. 

Using the Hamiltonian (2) the density matrix (1) can be represented as a function of 
parameters a = and (3 = 4I( ^flifk B T : 



P 



Z- 1 exp {-al z -(5U{9, <p) [si 2 z - P + \ {l\ - J 2 )] U + (9 ,<p)} . (7) 



The density matrix (7) can be employed to obtain information on the dependence of en- 
tanglement on the magnetic field, quadrupole coupling constant, orientation of the crystal 
principal axis in the laboratory frame, and temperature. 

Below we consider entanglement in a system of spins 3/2. A suitable system for studying 
by NQR technique is a high temperature superconductor Y Ba2Cu^Ojs containing the 63 Cii 
and 65 Cu nuclei with spin | possessing quadrupole moments Q = —0.211 x 10~ 24 cm 2 and 
—0.195 x 10~ 24 cm 2 , respectively [3^. There are two different locations of copper ions in 
this structure: the first are the copper ion sites at the center of an oxygen rhombus-like 
plane while the second one is five-coordinated by an apically elongated rhombic pyramid. 
The four-coordinated copper ion site, EFG is highly asymmetric (rj > 0.92) while the five- 
coordinated copper ion site, EFG is nearly axially symmetric (77 = 0.14). The quadrupole 
coupling constant (eQqzz) of 63 Cu in the four-coordinated copper ion site is 38.2 MHz and 



in the five-coordinated copper ion site is 62.8 MHz 



32|. 



III. REDUCED DENSITY MATRIX AND MEASURES OF ENTANGLEMENT 



An important measure of entanglement is the concurrence 33J. The concurrence C is 



usually used 



331 ] . For the maximally entangled states, the concurrence is C — 1, while for 



the separable states C = 0. The concurrence of a quantum system with the density matrix 
presented in the Hilbert space as a matrix 4 x 4 is expressed by the formula 331 ] : 



C = max ^ 0, 2v 



4 

1=1 



(8) 
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where v = max {ui, U2, ^3, ^4} and z/j (z = 1, 2, 3, 4) are the square roots of the eigenvalues of 
the product 

R = PredPred, (9) 

where p re d = Tr par u a i (p) is the reduced density matrix, Tr part i a i (■••) is the partial trace, and 



Pred 



(10) 



G 



(11) 



where p re d (a, (3) is the complex conjugation of the reduced density matrix p re d and 

/ -l\ 
10 
10 
\-l / 

An another important measure of entanglement is the entanglement entropy In order to 
implement this measure of entanglement we first map the Hilbert space of a spin 3/2 which 

bert space of two spins 1/2, Sa and S B (sSl- The entropy of 



is four dimensional onto the Hi 

2 



entanglement is defined as 



35|: 



E = -Tr (p A log 2 p A ) = -Tr (p B log 2 p B ) ■ 



(12) 



Here pa is the partial trace of pa = Tr B (p) over subsystem B, and p B has a similar meaning. 
This entropy is related to the concurrence C by the following equation 33] 



E (x) = —x log 2 x — (1 — x) log 2 (1 



(13) 



where x = |(l + ^l + C 2 ). 

The numerical simulations of entanglement of the spin states are performed for special 
case 1 = 3/2 using the software based on the Mathematica package. In the equilibrium the 
states of the system determined by Eq. (7) are separable without applying external magnetic 
filed (a = 0) at any temperature. At large temperature (/3 < 1) and low magnetic field 
strength (a < 1) the concurrence is zero. Entanglement appears in the course of temperature 
decrease and increasing the magnetic field. 

The concurrence and the entanglement of formation dependences on the magnetic field, 
quadrupole coupling constant, and temperature is qualitatively independent of the angles. 
As example, we present below the concurrence as a function of the parameters a and /3 



for the concurrence maximum for a spin in the five-coordinated copper ion site (77 = 0.14) 
at 9 = 0.94 and (p — (Fig. 1). The concurrence and the entanglement of formation 
increase with the magnetic field strength and inverse temperature and reache their maximum 
value. Then the concurrence and the entanglement of formation decrease with increasing 
the magnetic field strength (Fig. 2). Another dependence of the concurrence and the 
entanglement of formation on temperature is observed (Fig. 3). At a high temperature 
concurrence and the entanglement of formation are zero. With a decrease of temperature 
below a critical value the concurrence and the entanglement of formation monotonically 
increase till a limiting value. The critical temperature and limiting value are determined by 
a ratio of the Zeeman and quadrupole coupling energies, a/ (3. 



IV. DISCUSSION AND CONCLUSIONS 



The obtained results show that the entangled states can be generated between the states 
of a single nuclear spin |. From a point of of view of quantum information processing the 
considered system is isomorphic to a system consisting of two dipolar coupling spins ~. The 
same quantum logical gates can be realized using the both systems. Therefore the obtained 
entanglement can be considered as entanglement between qubits formed by states of a single 
particle. On the other hand, a single spin 3/2 isomorphic to a system consists of two dipolar 
coupling spins 1/2 and entanglement between the states of a spin 3/2 can be considered as 
entanglement between two effective spins 1/2. The behavior of entanglement of a spin 3/2 
is very similar to that for the system consisting of two dipolar coupling spins | [37]. It was 
obtained that in zero magnetic field the states of the both spin systems, two spins | and spin 
|, are in separable states. These systems become entangled with increasing the magnetic 
field. Then, with a further increase of the magnetic field the spin states of the both systems 
tend to a separable one. 

It has recently been shown that, in a system of nuclear spins s — 1/2, which is described 
by the idealized XY model and dipolar coupling spin system under the thermodynamic 



equ i 



[37|, 



ibrium conditions, entanglement appears at very low temperatures T w 0.5 -r- 0.3 /iaK 
In a non-equilibrium state of the spin system, realized by pulse radiofrequency 



irradiations, estimation of the critical temperature at which entanglement appears in a 
system of spins \ gives T < 0.027 K [39]. The calculation for 63 Cu in the five-coordinated 
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copper ion site of Y BCL2CU3O7-5 at a//3 = 1, rj — 0.14 and eQq zz = 62.8 MHz, gives that 
the concurrence appears at (3 = 0.6 (Fig. 4). This (3 value corresponds to temperature 
T w 5 niK. This estimated value of critical temperature is by three orders greater than the 
critical temperature estimated for the two dipolar coupling spins under the thermodynamic 
equilibrium [371 ] . 

In conclusion, performing investigation has shown that entanglement can be achieved by 
applying a magnetic field to a single spin 3/2 at low temperature. Concurrence and the 
entanglement of formation depend on the orientation between the external magnetic field 
and PAF axes. At 9 = and tt the states are separable at any conditions. At 77 = 0.14 the 
concurrence and the entanglement of formation reache their maximum value at 9 = 0.94 
and tp — and n . This effect can open a way to manipulate with the spin states by a 
rotation of the magnetic field or a sample. 
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Figure Captions 

Fig. 1 (Color online) The maximum concurrence as a function of the parameters a and 
(3 at r] = 0.14, 9 = 0.94, ip = 0. 

Fig. 2 (Color online) Concurrence (a) and the entanglement of formation (b) vs. magnetic 
field at T = const for various quadrupole interaction constants: black solid line - /3 — 2; red 
dashed line -j3 — 6 ; green dotted line -(3 = 8; blue dash-doted line - j3 — 12. 

Fig. 3 Concurrence (a) and the entanglement of formation (b) as a function of tempera- 
ture at | = 0.5 (black solid line), jf = 1 (red dashed line), and | = 2 (blue dotted line) at 
7/ = 0.14, 6> = 0.94, if — Temperature is given in units of 4J f 2 ^ z 1 f. . 
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